#include "NewtonInterp.H"
using namespace std;
int main(int argc, const char **argv)
{
  cout << "-----------------------Test NewtonIntp------------------------" << endl;
  NewtonInterp Solver1;
  RealBox b2 = {{0, 0}, {1, 1, 7, 42}, {2, 128, 448}};
  InterpConditions IntpConds2(b2);
  Solver1.ResetIntpConds(IntpConds2);
  cout << "//Test QuickInterp" << endl;
  Solver1.QuickInterp();
  Solver1.Display();
  cout << "//Test AddIntpSiteThenInter" << endl;
  Solver1.AddIntpSiteThenInterp({-1, -1, 7});
  Solver1.Display();
  cout << "//Test NevilleAitken" << endl;
  InterpConditions IntpConds3({{0, 0}, {1, 1}, {-1, 1}});
  cout << Solver1.NevilleAitken(-10, IntpConds3) << endl;
  return 0;
}